A Pair of Interacting Spins Coupled to an Environmental Sea : 
Dissipation and Mutual Coherence 

M. Dube and P.C.E. Stamp 
Physics Department, University of British Columbia, 6224 Agricultural rd., 
Vancouver, BC, Canada V6T 1Z1 



We consider the quantum dynamics of two spin-1/2 sys- 
tems, each coupled to a bath of oscillators, so that a bath- 
mediated coupling is generated between the spins. We find 
that the interaction destroys any coherent motion of the 2 
spins, even if the coupling of each spin to the bath is quite 
weak, unless the interaction is extremely small. This is be- 
cause the dynamic mutual bias between the spins blocks any 
coherent transitions between nearly degenerate states. In 
many quantum measurement operations this means that de- 
coherence effects will be much stronger during the actual mea- 
surement. 

PACS numbers: 



In this paper we study the theoretical problem of a 
Pair of Interacting Spins Coupled to an Environmental 
Sea (we use the acronym "P.I.S.C.E.S." for brevity). The 
2 spins will each be described as 2-level systems, and 
we shall concentrate on environments that can be rep- 
resented using a Feynman- Vernon oscillator bath [jl]]|] 
model - at the end of the paper we shall simply remark 
upon some analogous results obtained when the envi- 
ronment can be represented as a spin bath [||JJ]. Our 
PISCES problem is clearly a generalisation of previous 
models like the "spin-boson" model, describing a single 
spin coupled to a sea of oscillators - this model has been 
widely studied in connection with the Kondo problem 
H§, quantum diffusion [§, SQUID tunneling [|, and 
magnetic grains m, and extensively reviewed (see, eg., 
refs. |,0). 

The interest of the PISCES model is in the way it re- 
veals a dual role of the environment - whilst it causes 
decoherence and dissipation in the motion of each spin, 
one also expects that it might cause some mutual coher- 
ence between them. For very weak mutual coupling (as 
in the problem of two 2-lcvcl atoms coupled radiatively) 
this mutual coherence is fairly well understood. How- 
ever, as we shall see, a radical change occurs once the 
coupling exceeds a very small value, to be specified be- 
low. One then finds quantum relaxation behaviour of the 
2 spins, even at temperature T = 0. Since the majority 
of PISCES problems, in condensed matter and elsewhere, 
have couplings exceeding this value, our results have im- 
portant practical application. The PISCES model is also 
a primitive model for quantum measurements, in which 
both dissipation and the back-reaction of the measuring 
apparatus on the measured system are included. 

In cases where the environment is modelled by an os- 



cillator bath, we will start from an effective Hamiltonian 
Hpisc es = Hq + Hi n t , between spins located at position 
i"i and r2, where 



Ho = ™(Ai7f + A 2 r|) + 1 £ m k (±l + (1) 
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where rj* are the Pauli operators for each spin (J = 1,2), 
and {xk} are the bath oscillators. In this paper we are re- 
ally interested in the indirect bath-mediated interaction 
between T\ and T2, but we include also a direct interac- 
tion K a p(v) 1 with r = ri — r 2 ; K a p(r) can be considered 
as modelling the static part of some other field-mediated 
interaction [jll]. We assume that k is a momentum in- 
dex, describing delocalised bath modes; we also ignore 
non-diagonal couplings (ie., couplings in rf or 7\ y ), which 
have a less important effect on the physics ^,f7|. Clearly 
our model could be generalised, but at the expense of 
great complexity, and we believe that (Q) and @ bring 
out most of the essential physics. 

The interactions with the bath are characterised by 
spectral functions Jy (r, u>) whose diagonal elements 
Jj{oj) = Jjj(uj), referring to self-interactions of Tj via 
the bath, have the Caldeira-Leggett form [||: 



7T v-^ C 



C?)i 



-5{uj - uj k ) 



(3) 



The off-diagonal elements have a complicated form if 
|r| is large, because of retardation effects Q; these be- 
gin to be important once > Wk/|r|, where is the 
velocity of propagation of the k th oscillator mode. Here, 



we shall ignore such effects (but see ref. 13 1) which are 
most important in the weak-coupling limit, so that 
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where D\J0) is the static limit of the oscillator bath prop- 
agator |p2| . 

We shall consider 2 forms |L4| for Jy((r, w). The 
Ohmic form Jy(r, w) = cjrjij e~ u ^° , with r]j = rjjj a 
local friction coefficient and 7712 ~ (r^i^^V^r), is 
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appropriate to eg., electrons in Fermi liquid. The su- 
perohmic J#(r, uj) = g i3 (r)Q D (uo/Q D ) m e -"V B applies 
to phonons or photons, with m > 3 in 3 dimensions 
§; now g 12 ~ (gi3 2 ) 1/2 V (r). For electrons |V e (r)| - 
cos(2 k p r) / (kpr) 3 (the RKKY form) in 3 dimensions; for 
phonons |V^(r)| ~ (ao/r) 3 , where ao is the lattice spac- 
ing. The values of the upper cut-offs Oo and 0_d are 
typically Slo ~ eV for electrons, and ©£> ~ 10~ 2 eV" for 
phonons. 

We also define a mutual bias function e(r) as 



e(r) 
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with sign such that positive/negative e(r) corresponds 
to ferromagnetic/antiferromagnetic coupling between T\ 
and t 2 . Defining Qy(r) = rjij (r)qoiqoj /2irh, one then has 
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with e proportional to the upper cut-off frequency in both 
cases pif . 

A complete study of this problem calculates the 16 ele- 
ments of the 2-spin density matrix p (n, t 2 , t; t\, t 2 , t'). 
To do this we use standard influence functional tech- 
niques the total influence functional F{n, t 2 } can 
be factorised as F = Fix F 22 Fi 2 , where Fn and F 22 are 
identical to the single spin-boson functionals, and F\ 2 
contains all dynamic couplings via the bath. F\ 2 is a 
functional over the paths of both spins together, and 
is thus intrinsically more complex than F\\ and F 22 \ it 
cannot be evaluated in closed form for arbitrary values 
of the parameter in (S). The "dilute-blip" approxima- 
tion p| cannot be applied to F\ 2 in the very weak cou- 
pling limit because then one may have overlap in time 
between arbitrary configurations of each spin - thus this 
limit must be handled perturbatively Jig] , Away from 
this limit, the mutual bias e(r) plays an essential role. 
As a function of time the coupling energy in a path inte- 
gral contribution to pit) flips between ±e(r); once |e(r)| 
is large enough its effect is to suppress configurations in 
which both spins are simultaneously in blip states. One 
may then use a dilute-blip approximation for F\ 2 . The 
parameter regime in which this works depends on the 
precise form of Jy(r, w); we give the limits of validity 
below, for the ohmic and superohmic cases |]l5f 

There is no space here to give results for all matrix ele- 
ments of p (see ref. |1|); here we concentrate on P ++ {t) 1 
the probability to find the system in the state | tT) at 
time t, if at t = it was also in state | ff). Then for 
ferromagnetic coupling {a\ 2 > 0) we find 



P++(t) =A_+ A + e-^ +r ^ + B 4 
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(10) 



Ri 2 {T) = T\+Tl-'2T 1 T 2 [l - 2 tanh 2 (e/2fcT)] (11) 

with the rates T\ and T 2 being those for a spin-boson 
system in a bias e Q : 

Tj = A 2 / dt e- Q i W cos [Qtp (t) } cos [et/h] (12) 
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In this way the complete solution is seen to be one 
involving quantum relaxation of p(t) , with the relaxation 
rates (riTT2) and R± depending essentially on e/kT. To 
see how this works in practice, consider first the Ohmic 
case, where (2) 
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with the functions T(x) on the right-hand side being 
Gamma functions. If there were no interactions at all 
between t\ and t 2 , then at T = 0, the Tj would freeze 
for otj > 1, relax exponentially (with power law correc- 
tions) for 1/2 < ctj < 1, and show damped oscillations 
for a < 1/2; at any finite T, the Tj would relax expo- 
nentially for any value of the ctj (see ref. 0|). 

Switching on the interaction Ji 2 (y,uj) quickly changes 
this. Suppose, as above, that 012(1") is positive (ferro- 
magnetic). Then once e ^> fcT, ie., once «i2 ^> kT/Q a , 
so that tanh( e/2kT ) ~ —1, then P ++ (t) = 1 for all times 
of interests, ie., the 2 spins freeze each other ! More 
generally we see that at long times relaxes to 

j[l + tanh(ai2 Flo/kT)]. A similar freezing to the anti- 
ferromagnetic I tl) occurs if — ai2 3> kT/D,Q. Since f2o is 
usually large, this freezing can happen for very small cou- 
pling at low T. The results are valid if Aj <C |ai2|Ooj at 
any T, or for any ai2 if kT 3> Aj/ctj. In cases where the 
Aj describe tunneling matrix elements, as in SQUID's [|, 
quantum diffusion |Q] , or magnetic grains Q , the tunnel- 
ing rates are often rather small and so the results apply 
for very small couplings 012; the same will be true for 
many metallic glasses. 
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Similar results are obtained in the superohmic case, for 
which §§ 

Tj = (gojAj/e) 2 Jj{e) coth(e/2fcT) (16) 

In the common case where m = 3, this gives Tj ~ 
9oi5i2 r (3)A^/6D. If |gi 2 |e D > fcT, the spins will 
again be frozen; otherwise we get relaxation according to 
(0), provided Aj -C \gis \ ©d, for any T, or for any g±2 
if kT > Aj/^j e_D. Typically Aj < 6d and this relax- 
ation will be rather slow. This result shows that the bias 
e plays a crucial role - if there were no coupling between 
the spins, then they would not freeze even at T = 0, no 
matter how large the couplings g~\ and <? 2 . 

Some qualitative features of the results deserve special 
mention; in particular 

(i) It is quite common in phcnomcnological studies of 
this problem to simulate the environmental effects on 
the 2 coupled spins by stochastic terms |l6|]. In general 
such methods will give results quite unlike these (often 
with oscillatory behaviour in p(t) for large inter-spin cou- 
plings). The reason for the discrepancy is fundamentally 
that quantum influence functionals cannot be modeled 
by classical (c-number) noise source ]l7|] . 

(ii) The 3 relaxation times appearing in p(t) (see 
Eq.(n)) can De quite different, and their contributions will 
in general appear in different proportions in the 16 differ- 
ent elements of p(t). Fig.(|l|) illustrates this by plotting 
4 of these matrix elements (including P ++ (t)) calculated 
in the same way as (Q). 
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FIG. 1. Probabilities of occupation for a system with pa- 
rameters Ti = 0.5 sec" 1 , T 2 = 0.7 sec" 2 and l/2kT = -0.5. 
The system starts at t = in the state | ||); P a /3 is then the 
probability that at time t the system has rf = a and t% = f3. 

(iii) There is no trace left in our results of the "phase- 
transition" occuring at a = 1/2 and a = 1 in the spin- 
boson problem at T — 0, for the case of ohmic dissipa- 
tion M. The localisation transition at a = 1 is replaced 



here by a crossover to localised behaviour, at T = 0, as 
|e(r)| = 2f2 ai 2 (r) becomes greater than the larger of 
Ai and A 2 . We can understand this physically as "de- 
generacy blocking" , ie., the mutual bias blocks tunneling 
by destroying the degeneracy between the initial and fi- 
nal states. There is also a finite temperature crossover to 
localisation, for a finite coupling e(r) ^> Aj, as one low- 
ers T; this occurs for kT ~ 2 fio 012(1"), for ferromagnetic 
coupling, and kT ~ — 2 Slo 0:12(1"), for antiferromagnetic 
coupling. Analogous crossovers will also exist for super- 
ohmic coupling. 

A surprising feature of the results (at least for us) is 
that our solution (|7]) does not allow any mutually coher- 
ent oscillations of the two spins. Thus, eg., in the ohmic 
problem, if 01,02 < 1/2, so each spin exhibits damped 
oscillations if 012 = 0, one might have expected to find 
that if Ai ~ A2, switching on 012 could cause some 
kind of "frequency locking" between the 2 spins. Instead 
one simply finds that the mutual bias destroys the oscil- 
lations, once we are out of the very weak perturbative 
limit. A corollary of this is that if we start off with the 2 
coupled spins and no bath, then the effect of coupling in 
the bath is to destroy coherent oscillations of each spin 
far more quickly than one would expect from the known 
destructive effects of the bath on single spin coherence 

I- 

We have no space to discuss detailed applications 
of these results to particular systems (particularly as 
Ji 2 (r,(x>) is often long range in r, so that for a system 
containing many 2-level systems, it is not sufficient to 
simply examine the 2-spin correlation functions). How- 
ever our main result, that the coherence-destroying role 
of the bath is greatly exaggerated once mutual coupling 
between 2 (or more) systems is allowed, seems to be of 
general relevance. In particular it shows that during a 
quantum measurement, decoherence effects will be much 
stronger during the measurement operation. Consider, 
eg., a SQUID exhibiting damped coherent oscillations 
H for a < 1/2. It is clear that when this SQUID is 
coupled to a measuring apparatus (eg., another SQUID 
[jUI ) modelled also as a damped oscillatory 2-level sys- 
tem), the coherence will be destroyed unless the mutual 
coupling 012 is extremely small, ie., coherence will only 
persist whilst system and the apparatus are decoupled. 
A detailed analysis of this example [jl3| takes account of 
the coupling of each SQUID to the same set of EM field 
modes (whose detailed form will depend on sample ge- 
ometry) . The inductive part of the interaction which is 
responsible for the measurement is generated, as in the 
theory of this paper, by this coupling to the EM bath. 
One must also account for the resistive coupling of each 
SQUID's flux coordinate to separate quasiparticle baths. 

Finally we note that application of our results to many 
real systems will be complicated by the effects of the 
"spin-bath" of nuclear spins and paramagnetic impurities 
HQ, since these spins (a) add a further bias field acting 
on n and T2, and (b) can be flipped by a transition of 
r 2 or -r 2 , thereby causing further decoherence. The bias 
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field depends on the polarisation of the nuclear spins, 
and is particularly important if the system is magnetic. 
However it can even play a role in SQUID's, because 
although the SQUID flux couples only very weakly to 
the nuclear spins, it couples to a lot of them, and the 
total energy can easily exceed Aj by several orders of 
magnitudes. 
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